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PREFACE 


In  conjunction  with  RAND's  study  of  Defense  Against  Submarine- 
Launched  Ballistic  Missiles  for  the  Advanced  Research  Projects 
Agency,  background  investigations  of  theoretical  methods  for  calcu¬ 
lating  the  performance  of  nondirectional  passive  sonobuoys  are 
being  conducted.  This  particular  investigation  is  directed  toward 
detecting  broadband  signals  which  might  be  generated  by  certain 
types  of  submarines. 

The  author  is  an  Associate  Professor  in  the  Department  of 
Engineering  and  Applied  Science,  Yale  University,  and  is  a  consultant 
to  The  RAND  Corporation. 
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SUMMARY 

This  Memorandum  deals  with  some  of  the  problems  arising  in  the 
passive  detection  of  submarines  by  single,  nondirectional  hydrophones. 
The  signal  emitted  by  the  submarine  is  assumed  to  be  a  broadband 
noise  whose  characteristics  are  similar  to  those  of  the  background 
noise.  Hence,  detection  is  possible  primarily  because  of  the  increase 
in  noise  power  caused  by  the  presence  of  a  target.  The  effects  of  the 
motion  of  the  submarine  past  the  hydrophone  and  of  the  uncertainty  in 
the  background  noise  level  are  analyzed.  Although  the  observation 
time  for  a  slowly  moving  submarine  is  greater  than  that  for  a  fast 
submarine,  it  is  found  that  if  there  is  uncertainty  about  the  back¬ 
ground  noise  level,  detectability  is  sc-utially  independent  of  observa¬ 
tion  time.  Hence,  a  faster  and  therefore  noisier  submarine  is  more 
easily  detected  than  a  slower  and  quieter  one.  The  major  factor 
limiting  the  detection  range  is  the  noise  uncertainty;  but  if  the 
background  noise  is  stationary,  the  fact  that  the  signal  produced  by 
the  moving  target  is  a  nonstationary  noise  can  be  used  by  the  detector 
to  estimate  the  background  noise,  and  therefore  to  increase  the 
detectability. 
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I.  INTRODUCTION 

The  sound  emitted  by  a  submarine  is  assumed  to  consist  of  a 
number  of  line  frequencies  in  the  lower- frequency  part  of  the  signal 
band  and  a  broadband,  stochastic  type  of  signal  over  most  of  the 
higher- frequency  band.  In  this  Memorandum,  detection  based  on  the 
broadband  component  of  signal  is  considered.  For  certain  submarines 
the  amplitude  of  the  line  components  may  be  very  low  and  their  effect 
on  detectability  is  therefore  small.  Alsc ,  the  line  components  and 
broadband  components  are  probably  additive  in  their  effect  on 
detectability,  and  .he/  uan,  therefore,  be  considered  separately. 
Finally,  by  considering  only  the  broadband  spectrum,  a  scmewhat 
pessimistic  result  is  obtained  which  can  be  improved  by  processing 
the  line  components  in  an  optimum  fashion. 

This  analysis  is  preliminary  In  nature  and  its  purpose  is  to 
provide  rough  order-of-magnitude  estimates  of  detectability.  It  is 
based,  therefore,  on  a  number  of  simplifying  assumptions  concerning 
the  nature  of  the  target  signal,  background  noise,  and  transmission 
loss .  Specifically,  it  is  assumed  that  both  the  noise  and  the  signal 
measured  at  the  target  have  a  Gaussian  amplitude  distribution  with 
zero  mean  and  known  spectral  shape  (although  not  necessarily  known 
spectral  level)  and  that  they  are  stationary.  For  most  of  the 
analysis,  the  transmission  loss  is  assumed  to  be  inversely  proportional 
to  the  square  of  the  distance  between  source  and  receiver;  however, 
other  transmission- loss  curves  are  briefly  considered.  Any  effect  of 
the  transmission  medium  on  the  spectral  shape  of  the  signal  is  ignored. 
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The  target  submarine  is  assumed  to  travel  along  a  straight  line  at 
constant  velocity,  and  it  ..  assumed  that  this  velocity  is  slow 
enough  so  that  Doppler  shifts  are  of  negligible  importance,  at.  least 
to  the  broadband  detection  system. 

With  these  assumptions,  the  detection  problem  is  the  problem  of 
detecting  a  Gaussian  signal  in  a  Gaussian  noise  background.  This 
problem  has  been  studied  extensively,  (Ref.  1,  Chaps.  18-20,  and  Ref.  2, 
Chap.  XI)  and  the  pertinent  results  are  presented  in  Appendix  A. 
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II.  OPTIMUM  DETECTOR  FOR  DETECTING  A  MOVING  TARGET  IN  A  NOISE 
BACKGROUND  HAVING  KNOWN  STATISTICAL  PROBERTIES 

The  signal  received  by  the  sonobuoy  is  given  by 

x(t)  -  s(t)  +  n(t)  (1) 

where  s(t)  i3  the  signal  that  would  be  observed  if  there  were  no 
noise,  and  n(t)  is  the  noise.  In  the  present  case,  where  a  target 
submarine  is  assumed  to  travel  past  the  buoy,  the  signal  is  ampli¬ 
tude  modulated  and  is  therefore  in  the  form 

s(t)  -  £(t)  y(t)  (2) 

where  y(t)  is  a  stationary  (i.e.,  unmodulated)  stochastic  signal  and 
f(t)  is  the  deterministic  dimensionless  amplitude  modulation  resulting 
from  the  change  in  transmission  loss.  For  the  geometrical  situation 
shown  in  Fig.  1,  and  if  the  signal  power  is  inversely  proportional  to 
the  square  of  the  distance. 


Fig.  1— Target  and  sonobuoy  geometry 
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f  (t)  -  “=— * - (3) 

!i  +  (vt) 

where  k  is  a  constant  of  proportionality,  h  is  the  minimum  distance 
from  the  target  to  the  receiver,  and  v  is  the  target  velocity.  The 
time  t  ■  0  is  arbitrarily  chosen  as  the  time  at  which  the  target-buoy 
distance  is  a  minimum. 

A  nonstationary  signal  of  the  type  given  in  Eq.  (2)  is  most 
conveniently  expressed  in  terms  of  time  samples.  Thus,  the  received 
signal  x(t)  is  represented  by  an  n-dlmenslonal  vector 

il  "  fx(t1)»  x(t2),  -  x(tn)]  (4) 

where  the  prime  indicates  matrix  transposition.  According  to  the  sampling 
theorem,  if  the  signal  bandwidth  is  W  cps,  samples  taken  at  time 
intervals  separated  by  ^  seconds  represent  the  signal  completely  ex¬ 
cept  near  the  ends  of  the  observation  interval,  where  sampling  intro¬ 
duces  a  small  error.  For  an  observation  interval  of  T  seconds  the 
dimension,  n,  of  the  sample  vector  X  is  therefore  approximately  2  TW. 

The  optimum  detector  is  known  to  be  a  likelihood-ratio  detector 
which  decides  between  the  presence  or  absence  of  the  target  by  computing 
the  likelihood  ratio  corresponding  to  the  received  signal  and  compar¬ 
ing  it  with  a  preset  threshold.  The  decision  that  a  target  is 
present  is  made  if  the  threshold  is  exceeded.  The  threshold  is 
normally  adjusted  to  provide  a  specified  false-alarm  rate,  and  it 
can  therefore  be  regarded  as  a  function  of  the  false-alarm  rate. 

★ 

See  Ref.  1,  Section  4.2. 
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In  the  following  discussion  it  is  assumed  that  the  background 
noise  level  is  precisely  known.  It  is  shown  in  Appendix  A 
that  for  small  signai-to-noise  ratio  the  likelihood- ratio  receiver 
is  equivalent  to  a  device  that  computes  the  test  statistic 

u  121  X  (5) 


where  2  and  P  are  covariance  matrices  of  &  corresponding  to  noise 
only  present  and  to  signal  only  present  respectively,  i.e. 


2“ 


N 


(6) 


and 


The  symbol  ^  ^  represents  the  statistical  average  of  the  quantity 
within  the  bracket  conditional  on  the  hypothesis  indie  ted  by  the 
subscript.  The  statistic  U  is  ceopared  to  a  threshold  Uo-  For  low 
signai-to-noise  ratir,  it  is  shown  in  Appendix  A  that  the  figure  of 
merit  of  the  optimum  detection  system  is  the  quantity 

d  •  v4  tr[(P  Q'1)2]  <*) 

where  tr(  )  denotes  the  trace  of  the  matrix. 

In  order  to  consider  the  effec-:  of  the  distance  modulation, 
suppose  for  the  moment  that  both  the  spectra  of  the  signal  y(t) 
and  of  the  noise  background  are  fiat  ("white")  over  the  frequency 
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band  0  £  f  £  W  and  that  they  vanish  for  other  frequencies.  In 
practice  it  is  found  that  both  signal  and  noise  spectra  fall 

(3) 

off  at  approximately  the  second  power  of  frequency.  However,  it  can 
be  shown  that  the  results  are  not  materially  affected  by  the  exact 
spectral  shapes  as  long  as  they  are  approximately  the  same  for  signal 
and  noise,  and  as  long  as  the  modulation  is  very  slow  relative  to 
^  ,  as  would  be  the  case  here.  If  white  noise  is  assumed,  the 
individual  samples  are  uncorrelated,  and  the  covariance  matrix  £  of 
the  noise  background  is  simply 


a-Ni  (9) 

2 

where  N  is  the  average  noise  power  (in  units  such  as  nbar  )  and  I 
is  the  n- dimensional  unit  matrix.  Similarly,  the  covariance  matrix 
of  the  sample  vector  Y  corresponding  to  the  unmodulated  stochastic 
signal  y(t)  defined  in  Eq.  (2)  has  the  form 


(10) 


where  S  is  the  average  signal  power  prior  to  modulation.  As  a  result 
of  the  modulation,  the  covariance  matrix  P  has  the  form 
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Hence,  the  figure  of  merit  d  of  Eq.  (8)  becomes 


d 


which,  by  substitution  of  Eq.  (3),  becomes 


(12) 


d 


1  n 

r 

IV 

1 

2  L 

,2  „  .2 

V  i-1 

h  +  (vt£) 

(13) 


The  summation  is  most  easily  evaluated  by  converting  it  to  an 
integral,  which  is  permissible  if  v2(ti  -  t£  L)2/h2  «  1  as  would 
normally  be  true  here.  Then 
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i-1 


~i 


2  2 
h  +  (v  t.) 

l/ 


T/2 

x  r  . 

dt 

T/2 

M  Of  7  /  - 

dt 
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2  2 
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|  J 
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where  At  -  t£  -  t^  »  —  by  using  the  sampling  theorem.  The  integration 
is  straightforward  and  yields 


d  = 


vT 

2h 


vT 

2h 


+  tan 


-1  vT 
2h 


(14) 


It  is  clear  that  d  is  maximized  by  letting  T 


00 ;  in  this  case 
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d 

max 


1.253 


(15) 


However,  If  T  is  only  2h/v,  which  means  that  the  observation  time  is 
defined  for  the  angle  9  in  Fig.  1  to  vary  from  -45°  to  +45°,  then  d 
has  reached  98  per  cent  of  its  maximum  value.  Thus  T  ■  2h/v  might 
be  used  as  a  practical  definition  of  the  time  that  the  target  is  within 
sonar  range. 

The  implementation  of  the  physical  device  that  forms  the  statis¬ 
tic  U  given  by  Eq.  (5)  fellows  directly  by  substituting  Eqs.  (9)  and 
(11)  into  Eq.  (5).  This  gives  the  result 

°  ■  i  i  f2(tt}  *i2 

2N  i-i 


ws 

N2 


f2(t)  [*(t>] 


2 


dt 


(16) 


where  '.he  conversion  of  the  sum  to  the  Integral  follows  in  the  same 
way  as  in  Eq.  (13).  Tills  result  implies  that  the  detector  consists 
of  a  square-law  device  followed  by  a  cross  correlator  or  filter 
matched  to  the  modulation  envelope  as  shown  in  Fig.  2.  It  is  clear 
that  the  system  cannot  be  constructed  unless  f(t)  is  known,  which 
means  that  h,  v,  and  the  time  at  which  the  target  is  closest  to  the 
buoy  be  known  a  priori,  i.e.,  before  the  target  appears.  Thl*  is 
obviously  an  unreasonable  requirement,  and  we  consider,  therefore, 
a  simpler  suboptimum  scheme. 


Fig. 2— Likelihood  ratio  detector 
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III.  SIMPLE  SUBOPTIMUM  DETECTOR  FOR  DETECTING  A  MOVING 
"  TARGET  IN  NOISE  WITH  KNOWN  STATISTICS 


The  suboptimum  detector  considered  here  consists  simply  of  an 
integrator  which  integrates  the  squared  received  signal  for  T  seconds, 
i.e.,  it  puts  out  a  test  signal 


T/2 


An  optimum  value  of  T  will  be  derived,  and  the  relative  magnitude  of  the 
detection  index  d  produced  by  this  system  will  be  obtained. 

As  shown  in  Appendix  A,  for  small  signal- to-noi  e  ratio  the 
detection  index  d  can  be  written  in  the  approximate  form 


d 


(18) 


where  and  are  respectively  the  mean  values  of  U  under  the 
hypothesis  that  signal  is  present  and  that  it  is  absent,  and  aQ  is 
the  rras  fluctuation  of  U  for  signal  absent.  For  the  suboptimum 
system  we  define  the  detection  index  in  the  same  way  with  the  m-'s  and 
a  corresponding  to  U  . 

If  signal  is  present,  X(t)  =  f(t)  y(t)  +  n(t),  and  therefore 
from  Eq.  (17) 


U 


y2(t)  +  2f(t)  y(t)  n(t)  +  n2(t)]  dt 


(19) 
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The  mean  value  n*  is  obtained  by  averaging  each  term  in  the  integrand 
separately.  The  cross-product  term  vanishes  because  y(t)  and  n(t) 
are  independent.  Then,  using  the  fact  that  y(t)  and  n(t)  are 
stationary  and  that  the  signal  and  noise  powers  are  S  and  N 
respectively 


* 


T/2 


1  f2(t)  dt 

+ 

—  m 

W  T  “ 
w  i  N 

-T/2  J 

(20) 


The  integration  of  f^(t)  given  by  Eq.  (13)  is  straightforward  and 
yields 


* 


2  2 
2k_  WSl 


hv 


N 


“'l  a  4  WT5 


(21) 


If  signal  is  absent,  f(t)  in  Eq.  (19)  is  zero,  and  therefore  the  first  ten 
in  5q.  (21)  vanishes;  hence 


* 


H  I  - 

N 


(22) 


Also,  it  is  easily  shown  that  for  Gaussian  noise  and  sample  size  n  *  2TV 


WTt 

N 


(23) 


Thus,  the  suboptimum  detection  index  d/  is 


*  *  _ 

-  14 1  *  ^o  kZ  S  J~  2h  .  -1 

■  - - r -  ■  —  —  VT  W  —  tan 

a  .  2  N  vT 

o  h 


vT 

2h 


(24) 
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This  has  been  plotted  in  Fig.  3,  together  with  the  optimum  detect  Lon 
index  given  in  Eqr  (14).  Figure  3  shows  that  the  optimum  value  of 
T  is  approximately  3  ~  .  With  this  value  of  T  the  suboptimum  detec¬ 
tion  index  reaches  about  90  per  cent  of  the  absolute  optimum  given  by 
Eq.  (15).  Figure  3  also  shows  that  the  index  is  not  very  sensitive  to 

relatively  large  changes  of  T  away  from  the  optimum  value;  thus  for  T 
h 

between  3/4  —  and  15  —  the  index  exceeds  60  per  cent  of  the  absolute 
v  v 

optimum. 

This  result  indicates  that  if  the  receiver  simply  integrates 
the  squared  received  signal,  the  integration  time  should  be  optimized 
relative  to  the  most  distant  target  that  the  system  can  reasonably 
be  expected  to  detect.  Then  for  a  target  at  a  smaller  range,  T  will 
be  too  large  to  be  optimum;  but  since  the  target  is  closer,  it  will 
be  more  easily  detectable  so  that  the  small  loss  in  optimality  of  T 
is  of  no  consequence. 


imum  detector 
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g. 3— Detection  indices  for  optimum  and  suboptimum  detectors 
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IV.  EFFECT  OF  PRACTICAL  SIGNAL  AND  NOISE  SPECTRA 


In  Sections  II  and  III  the  spectra  of  signals  and  noise  were 
assumed  to  be  white  in  order  to  simplify  the  analysis  of  the  problem 
of  motion  of  the  target  submarine  past  the  sonobuoy.  Since  the 
analysis  indicated  that  the  effect  of  the  distance  modulation  can  be 
approximated  quite  well  in  the  detector  simply  by  integrating  the 
received  signal  over  an  appropriate  length  of  time,  it  seems  reason¬ 
able  to  make  the  further  approximation  that  the  signal  received  from 
the  target  has  constant  power  during  the  observation  interval.  This 
approximation  makes  it  possib1:;  to  treat  both  signal  and  noise  as 
stationary  processes,  so  that  other  than  white  spectra  can  be  easily 
dealt  with. 

When  it  is  desired  to  consider  the  detection  of  stationary 

signals  with  complicated  spectral  properties,  it  is  most  convenient 

(4) 

to  expand  the  received  signal  in  a  Fourier  series 

m 

cos  tu^t  +  b^  sin  iw^t) 
k»l 


2tt 

where  **  Aud  »  —  with  T  the  observation  interval.  The  a^'s  and  b^ 's 

are  Gaussian  random  variables  with  zero  mean  value;  and  for  T  very  much 

larger  than  the  inverse  bandwidth  of  the  received  signal,  they  are 

approximately  mutually  independent,  i.e.,  a.b  »  0  for  all  k  and  l,  and 

k  l 


a  a  =  b.b  =0  for  k  ^  Their  mean-square  value  is  given  by 
k  }  k  l 


(4) 


G(fk)  Af 

\  G(V 


(26) 
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where  G  (f)  is  the  power  spectral  density  of  x(t)  and  f  ■  — .  If 
x(t)  consists  of  noise  only 

G  (f)  =  N  N  (f)  (27) 

o 

where  Nq  is  the  low-frequency  spectral  level  of  the  noise  which  is 
assumed  to  be  finite.  Thus 

Lim  N  (f)  «  1  (28) 

f~0 

If  x(t)  consists  of  signal  and  noise,  then  because  of  the  independence 
of  signal  and  noise 

G(f)  -  Sq  S(f)  +  M  N(f)  (29) 


where  S  is  the  low-frequency  spectral  level  at  the  receiver  (again 
o 

assumed  finite)  and  S(f)  is  the  normalized  spectral  density  of  the 


signal  with 


Lim 

f~o 


S(f)  -  1 


(30) 


The  reasons  for  normalizing  the  spectral  densities  are  discussed  in 
Appendix  B;  they  simplify  the  discussion  of  signal  with  unknown  power 
level. 


The  standard  theory  of  detection  reviewed  in  Appendix  A  can  be 

applied  to  the  formulation  in  this  section  by  considering  the  Fourier 

coefficients  a^  and  as  the  elements  of  the  sample  vector  X,  t.e.  , 

X  ■  fa, i  b, ,  a„.  b.  ....  a  ,  b  1  where  n  is  large  enough  so  that  all 
-  u  1  1  2  2  nn 

significant  frequencies  are  included.  Considering  Eqs.  (27)  and  (29), 
and  because  of  the  independence  of  the  coefficients,  the  matrices  P 


and  ^  become 
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l  -  So£  "  T 


S^) 


S(f  x) 


0 


s(f2) 


s(f2) 


0 


•S(fn) 


N 


2.  “  N  q  - 


T 


N(fx) 


N(f 


0 


N(f2) 


N(f2) 


0 


•N(fn) 


(31) 


(32) 


Then  the  normalized  test  statistic  considered  in  Appendix  B  (Eq.  (B-3)) 
becomes 


u 


2-3- 


-1 


El'1  £ 


T 

2 


n 


I  +  b*> 


s<fk> 

[N(fk) j2 


(33) 


It  is  easily  seen  that  u  is  essentially  the  output  of  a  system  of  the 
type  shown  in  Fig.  4,  where 

|H  (u))l2  =  (34) 

tN  (a,)  r 

is  identical  to  an  optimum  filter  discussed  by  Zadeh  and  Ragazzini.^ 


Fig. 4— Optimum  detector 

This  can  be  shown  by  computing  the  output  of  the  squarer  which  is  given 
by 
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v(c)  -  fe  afc  |HUk)|  cos  (o)kt  +  6k)  *  ^k  |H<uk)|  si"  +  V}j 
where  0.  Is  the  phase  angle  of  H(u^)  .  Equation  (33)  then  results  since 
for  large  T  the  Integrator  output  is  Tvdc  where  vdc  is  the  dc  component 
of  v(t).  H(uj)  must  be  stable  and  physically  realizable,  but  this  re¬ 
quirement  causes  no  difficulty  in  practice. 

Expressions  for  the  false-alarm  probability  and  the  probability 
of  true  detection  have  been  derived  in  Appendix  B  for  the  case  of  unknown 
but  small  signal  power  and  are  rek  roduced  here  for  convenience 
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where  passage  to  the  integral  is,  as  usual,  permi.sible  if  T  is  very 
large  compared  to  the  reciprocal  of  the  bandwidth. 

Example  broadband  spectra  for  submarine  targets  and  bacKground 
noise  are  shown  in  Fig.  5.  The  spectra  obtained  from  actual  measure¬ 
ments  may  differ  considerably  from  these,  and  their  exact  shape 
depends  on  a  number  of  factors  in  addition  to  target  velocity  and  sea 
state. ^  However,  Eqs.  (35)  to  (38)  indicate  that  small  differences 
in  the  shapes  of  signal  and  noise  spectra  have  negligible  effect  on 
either  the  false-alarm  or  detection  probability.  The  normalized 
spectra  can  therefore  be  approximated  by 


S(lu) 


2 

ml 

2  2 
u>  +  uu  ^ 


(39) 


N(««) 


<JU, 


2. 

+  u>2) 

27 

+  »0> 


(40) 


Both  signal  and  noise  spectra  fall  off  with  the  second  power  of 
frequencies  above  a  frequency  or  respectively,  and  the  noise 
spectrum  levels  off  at  a  high  frequency  uj^  because  of  the  presence 
of  locally  generated  white  noise. 

The  integrals  of  Eqs.  (37)  and  (38),  after  substitution  of 
Eqs.  ^39)  and  (40),  are  in  a  standard  form  that  is  tabulated,  for 
example,  in  Ref.  7.  The  result  is 


T  tM1  0>2  (0)1  tt)2  +  0io) 

4  2  ,  .  N 

%  °'-i  +  V 


(41) 
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IP 


Fig. 5— Example  spectral  densities  for  target  and  for  noise 
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and 


T“l  “2  ,  3 

4  4  '-U)l 

4tu  Co)1  +  ui„) 

0  1  /. 


(u>2  +  %) 


2  4  2  2  2  2 

+  4u),  tu,  tu  +  2u),  w,  uj  (w1  +  2u)  ) 

1  L  o  1  L  o  1  o 

+  0^2  (to*  +  u^)2]  (42) 

Typical  values  for  ,  uj,  .  and  are 

uj  ■  2rr  x  50 
o 

uj  -  2tt  x  1000  (43) 

1 

oj2  -  2tt  x  20,000 

Hence  m0  >  >  tu,  >  0)  ,  and  Eqs.  (41)  and  (42)  can  be  approximated  by 

2  1  O 


tr  (£a  S 


T_ 

2n 


/  -U2 
tr  (££  ) 


!_ 

2tt 


(44) 


(45) 


The  false- alarm  and  detection  probabilities  can  now  be  evaluated  by 
substituting  Eqs.  (44)  and  (4i.v  .n  Eq^.  (35)  and  (36).  It  must  be 
noted,  however,  that  the  low-frequency  spectral  level  Sq  of  the 
signal  refers  to  the  level  at  the  receiving  sonobuoy,  whereas  the 
level  given  in  Fig.  5  is  referenced  to  a  distance  of  one  yard  from  the 


target 
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If  the  transmission  loss  is  inversely  proportional  to  the 
square  of  the  distance,  then  by  Eq.  (3) 


S 

o 


,  2  ,  s2 
h  +  (vt) 


where  is  the  low-frequency  spectral  level  of  the  target  relative 
to  one  yard  and  h  is  the  minimum  range  in  yards.  The  assumption 
that  the  target  signal  strength  at  the  sonobuoy  is  constant  during 
the  observation  time  implies  that  vt  «  h  so  that 


S  ss  S.  - r 
0  1  h2 


06) 


Thus  the  false-alarm  and  detection  probabilities  become 


i  -  e 


1 

2 

1 

2 


;0(z) 

2  <y 


2  8  <  V 


(47) 

(48) 


where 


2  u 


a 


n  riu. 

o  2 


-a  - 

V 


(49) 


and  where 


1  - -  /  \  (  k  i  I  ^1 

r  Mrs)  4 

'  o/  \h  /  \'JbQ 


(50) 


2  2  -t2 

©  (Z)  -  7=-  f  e  dt 

V11 


* 
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The  magnitude  of  T  in  the  above  expression  has  been  considered 

in  Section  II  and  was  shown  to  depend  on  the  range  h.  The  detectable 

range  might  be  defined  as  the  value  of  h  that  makes  Zg  -  0  subject 

to  a  particular  false-alarm  rale.  For  fixed  false-alarm  rate 

-5 

is  a  constant,  k^;  specifically  for  a  -  10  ,  k^  •  3.03.  Letting 

T  ■  3h/v  according  to  the  results  of  Section  II,  we  can  solve  Eq,  (50) 

for  h  ,  i.e. 
max 


or 


(51) 


For  v  -  20  knots,  we  find  from  Fig.  5  that  S^  -  27  db.  Also, 

assuming  that  the  sea  state  is  1,  N  "  -24  db.  u>o,  u.’^,  and  ar®  given 

in  Eq.  (43).  Substituting  all  of  these  values  into  Eq.  (51)  with 

k  -3.03  results  in  h  of  about  600  mi,  and  T  is  on  the  order  of 
ot  max 

100  hr.  These  results  are  clearly  unrealistic,  and  they  indicate 
that  an  important  factor  has  been  omitted  from  the  analysis. 
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V.  EFFECT  OF  UNKNOWN  NOISE  LEVEL 

The  unreasonable  result  obtained  in  Section  IV  can  be  shown  to 
result  from  the  assumption  that  the  noise  level  is  known  precisely. 

In  order  to  see  this,  consider  what  actually  happens  in  the  detection 
system.  The  detector  is  shown  in  Fig.  4,  and  if  noise  only  is  present 
at  the  input,  it  is  clear  that  u  will  be  a  somewhat  random  ramp 
function  as  shown  ii.  Fig.  6.  If  signal  is  also  present,  then  the 
slope  of  the  ramp  is  slightly  greater,  as  shown.  Thus,  at  any  time 
T  it  is  possible  to  establish  a  threshold  which  will  usually  reject 


Fiy, 6— Detector  output  versus  time 


the  signal  hypothesis  when  it  is  incorrect  and  accept  it  when  it  is 
correct  (the  word  "usually"  is  used  in  a  statistical  sense).  It  can 
be  seen  from  Fig.  6  that  as  long  as  the  input  power  level  with  signal 
plus  noise  is  even  very  slightly  greater  than  it  is  with  noise  alone, 
the  average  slopes  of  the  two  rairns  are  different  and  they  continue 
to  separate  as  the  observation  time  is  increased. 

This  discussion  is  valid  if  one  can  assume  that  the  noise  power 
is  exactly  known.  Suppose,  however,  that  there  is  some  uncertainty 
about  the  noise  power  level.  Then,  in  order  to  achieve  a  specified 
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false-alarm  probability,  the  threshold  would  have  to  be  set  at  a 
high  enough  level  so  that  the  maximum  expected  noise  level  would 
yield  the  specified  false-alarm  rate.  If  then  the  noise  level  is 
actually  less  than  the  maximum  value,  u  would  usually  be  below 
the  threshold  for  both  signal  present  and  absent  unless  the  signal 
power  exceeds  the  uncertainty  in  the  noise  power. 

In  view  of  this  discussion,  one  can,  for  a  given  observation 
time  X,  define  a  minimax  range  h^  as  the  range  that  causes 

to  vanish,  with  N  taking  on  its  smallest  value,  for  a  Z  taking 

o  <y 

on  the  value  corresponding  to  the  specified  false-alarm  probability 

with  the  largest  value  of  N  . 

o 

Quantitatively,  suppose  that  the  noise  level  is  somewhere 
between  N  4-  AN  and  N  -  AN  .  Then  if  the  value  of  Z  of  Eq.  (49) 

OOOO  Ct 

is  fixed  at  K 

Of 


K 

a 


2u 

_ o _ 

(Nq  +  ANo)v/TS£ 


02) 


The  value  of  Z  to  be  used  in  Eq.  (50)  is  th:n 

a 


Z 

a 


(K  + 

a 


|yTu)2) 


N  4  AN  \ 
o _ o  \ 

N  -  AN 
o  o  ) 


(53) 


and  h  is  obtained  from 
m 
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A  set  of  curves  of  h  versus  AN  /N  is  given  in  Fig.  7  for  various  sea 

in  o  o 

Ttates,  and  for  a  target  velocicy  of  20  knots.  The  example  spectra 

of  Fig.  5  have  been  used  in  this  cc -'putation.  Values  of  uj  .  uj.  ,  and 

o  1 

11)2  are  those  givun  in  Eq.  (A3).  For  these  values,  and  for  AN  /Nq  >  .01, 

the  first  term  in  Eq.  (5A)  is  negligible,  and  a  very  good  approximate 

solution  for  h  is 
m 


h 

m 

k 


(55) 


The  optimum  time  of  observation  is  given  by  T  ■  which  for  a  20- 

knot  target  velocity  is  T  *=  .15  h  hr  if  h  is  expressed  in  n  mi.  Thus 

m  m 

the  optimum  time  can  be  read  from  Fig.  7.  According  to  Eq.  (55)  dif¬ 
ferent  target  velocities  affect  the  range  only  through  the  signal 
spectral  level  S^.  According  to  Fig.  5,  for  a  5-knot  target  ie 
about  10  db  below  the  value  for  a  20-knot  target.  Hen'-a,  h  is  .316 

IQ 

as  large  for  a  5-knot  target  as  for  the  20-knot  target,  «.nd  the  optimum 

value  of  T  becomes  ,6h  hours. 

m 

By  setting  N  +  AN  *  N  and  N  -  AN  -  N  .  ,  Eq.  (55)  can  be 
o  o  max  o  o  mm 

put  into  the  form 


h 

m 

k 


1 


N 


max 

I  , 
min 


-  1 
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Although  this  expression  is  equivalent  to  Eq.  (65),  it  is  more  con* 

venient  when  N  /N  .  is  large.  It  has  been  used  for  the  curves 
max  min 

shown  in  Fig.  8,  which  are  also  based  on  the  example  spectra  of 

Fig.  5  and  on  the  values  of  oj  ,  U).  ,  and  u;  from  Eq.  (43). 

o  i  i 

It  should  be  noted  that  if  the  approximate  expression  Eq.  (55) 

holds,  then  h  is  independent  of  T  since  T  enters  Eq.  (54)  through 
m 

the  square-root  factors  that  have  been  cancelled  out.  Thus,  although 
it  is  desirable  that  the  adjustment  of  T  be  approximately  correct  in 
order  to  optimize  the  detection  index  as  discussed  in  Section  III, 
the  fact  that  different  trrget  velocities  require  different  values  of 
T  would  not  have  to  be  considered  in  the  design  of  the  detection 
system.  T  could  be  determined  by  measurements  of  the  sea  state  (or 
noise  background)  if  the  expected  uncertainty  in  the  noise  measure¬ 
ment  is  known. 

In  practice,  it  appears  that  accuracies  of  i  1  db  are  about  the 

it 

best  that  can  be  expected.  These  accuracies  correspond  to  =  .1, 

and  therefore,  according  to  Fig.  7,  for  sea  state  0  the  raiuimax 
detection  range  for  a  20-knot  submarine  is  about  25  mi,  with  more 
normal  sea  states  yielding  ranges  of  about  8  mi  or  less. 

The  probability  of  detection  as  a  function  of  range  depends  on 
the  probability  distribution  of  the  actual  noise  level  within  its 
permissible  range.  The  only  case  considered  in  this  section  is  that 
this  distribution  is  uniform,  i.e.,  if  the  actual  noise  level  is  N^ 

★ 

Private  communication  with  J.  Kingsbury,  Navy  Underwater  Sound 
Laboratories,  New  London,  Connecticut. 
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Fig. 7— Minimax  range  versus  uncertainty  in  noise  ievel 


Nmo*  =  Maximum  zero  frequency  noise  level 
Nmin  =  Minimum  zero  frequency  noise  level 
Target  velocity:  20  kt 

Transmission  loss:  6  db  per  distance  doubled 


max 


Fig. 8 — Minimax  range  versus  noise  uncertainty 
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dN 

P  t«l>dNl  ■  -  No  '  aNo  <  N1  <  No  +  iNo 

o 


p  (N.,)  =  0  otherwise 


(56) 


The  computation  for  the  detection  probability  is  simplified  by  the 
fact  that  for  a  given  value  of  ,  (1-p)  decreases  very  rapidly  for 
small  changes  in  h  near  the  nominal  value.  If  T  is  one  hour  and  0)2 

g 

is  2tt  x  20,000  rad/sec,  Tu^  *  4.5  x  10  rad.  For  this  value  of  Tu^ 
it  can  be  shown  by  direct  computation  of  (1-0)  according  to  Eq.  (48) 
that  if  (1-0)  =  \  for  h  =  hQ,  (1-0)  1  for  h  -  .999  hQ  and  (1-0)  c*  0 

for  h  =  1.001  hQ.  Hence,  if  we  use  the  notation  p(D/Np  instead  of 
(1-0)  to  indicate  that  the  detection  probability  is  conditional  on  N^, 
we  have  approximately 


p  (D/Nx)  -  1  -  u  [h  -  h  (N,)]  (57) 

where  u  [  ]  is  the  unit  step  function  and  where,  by  an  argument 
similar  to  that  leading  to  Eq.  (55) 

h  (N  )  uk  ys7 

k  "  “  /f  +  AN"  N,  (58) 

o  v  o  o  1 


The  joint  probability  p  (D,  N^)  obtained  by  multiplying  Eqs.  (56) 
and  (57)  may  now  be  integrated  over  to  obtain  the  desired  marginal 
probability  density  P  (D).  The  integration  is  straightforward  ar.d 
yields 


P  (D) 


2  2 

uuj  Sl  YT 

2~2 

2 u)  h  AN 
o  o 


for 


(59) 


P  (D) 


for 


This  expression  is  plotted  in  Fig.  9.  It  is  see  i  that  the  probability 
of  detection  is  unity  for  ranges  less  than  the  minimax  value  and  that 
it  drops  off  rather  rapidly  for  larger  ranges. 

The  detection  range  h^  is  conznonly  defined  as  the  value  of  the 
range  for  which  the  detection  probability  p(D)  ■  if.  It  is  clear  from 
Eq.  (59)  that  this  is  given  by 


This  range  can  easily  be  read  from  Figs.  7  and  8  by  a  change  in  the 
ordinate  scale. 

It  might  be  noted  that  a  change  in  the  false-alarm  rate  has 
essentially  no  effect  whatever  on  the  minimax  or  detection  ranges, 
or  on  the  detection  probability.  A  change  in  the  false-alarm  rate 
results  in  a  change  in  K^,  but  the  term  in  Eq.  (54)  containing  K 
was  neglected  in  obtaining  Eq.  (55),  and  in  all  subsequent  equations. 
Hence,  the  results  of  Figs.  7,  8,  and  9  do  not  depend  on  the  false- 
alarm  probability. 

The  results  of  this  section  ere  easily  extended  to  transmission- 
loss  curves  other  than  the  simple  6  db  per  distance  doubled  curve 
considered  thus  far.  Suppose  that  Eq.  (46)  is  replaced  by  the  more 
general  equation 


S 
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si 


(61) 


Then,  if  in  Eq.  (54)  the  term  involving  K  is 

a 


again  ignored,  Eq.  (55) 


becomes 
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h  = 
m 


(62) 


Also,  Eq.  (59)  becomes 
P(D)  *  1 


i’(D) 


f  h,  ,  AN° 
for  2  2  ~ 


t  \  + 

Uj  \ 

-2.  i 


uo 
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(63) 


v,  AN 
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under  certain  conditions  the  transmission-loss  curve  appears  to  have 
approximately  an  8  db  per  distance  doubled  slope.  Then  g(h/k)  =  (b/k) 
The  minimax  range  corresponding  to  this  form  of  g(h/k)  is  plotted  in 
Figs.  10  and  11.  The  larger  transmission  loss  results  in  a  con¬ 
siderable  reduction  in  range,  as  might  be  expected.  As  before,  the 

detection  range  can  be  defined  by  h.  «* y/2  h  . 

d  m 

As  a  result  of  multiple  reflections,  g(h/k)  frequently  has  a  more 

complicated  form,  of  the  type  shown  in  Fig.  12,  where  the  peak  after 

the  first  dip  typically  occurs  at  10  to  20  n  mi.  It  is  clear  from 

Eq.  (63)  that  under  these  conditions  the  form  of  the  detection 

probability  curve  will  be  as  shown  in  Fig.  12  and  that  there  may  be 

several  widely  different  values  of  h  and  h,. 
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Fig.  10— Minimax  range  versus  uncertainty  in  noise  level  for 
8  db  per  distance  doubled  transmission  loss 
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Fig.  11 — Minimax  range  versus  noise  uncertainty  for 
8  db  per  distance  doubled  transmission  loss 
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VI.  UNCERTAINTY  OF  NOISE  LEVEL  DESCRIBABLE  BY  A  GAUSSIAN 
PROBABILITY  DENSITY  FUNCTION 


In  Section  V,  it  was  assumed  that  definite  upper  and  lower 
bounds  coula  be  placed  on  the  possible  value  of  the  noise  level,  and 
in  the  computation  for  the  detection  probability,  it  was  assumed 
that  all  noise  levels  inside  these  limits  were  equally  likely. 

Although  it  is  conceivable  that  certain  techniques  for  determining 
the  noise  level  might  result  in  noise  uncertainties  of  this  type,  it 
seems  more  likely  in  general  that  the  noise-level  uncertainty  would 
be  more  accurately  described  by  a  probability  density  function  such 
as  the  Gaussian  density  function,  for  .;hich  it  is  impossible  to  define 
upper  and  lower  bounds  on  the  noise  level.  Therefore,  the  concept 
of  a  minimax  range  as  defined  in  Section  V  is  not  applicable.  However, 
if  the  probability  density  function  of  the  noise  level  is  known,  it 
is  still  possible  to  determine  a  threshold  resulting  in  a  specified 
average  false-alarm  rate. 

If  the  probability  density  function  of  the  noise  level  is  p(N^), 
then  th  average  false-alarm  and  detection  probabilities  are, 
respectively 

CD 

a  -  J  a(V  PvN^  dNx  (64) 

-00 

00 

(1-B)  -  1  -  J  B(N1)  pCNj)  dN]  (65) 

“00 


where  ar(N^)  and  @(Np  are  the  conditional  false-alarm  and  miss  pro¬ 
babilities  under  the  assumption  that  the  noise  level  is  N^,  i.e., 
they  are  given  by  Eqs.  (47)  and  (48)  with  replaced  by  N^, 
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Oue  way  of  estimating  the  noise  level  is  to  make  a  long-time 
measurement  of  the  received  signal  power  and  to  use  the  average 
power  as  the  estimate  of  noise  power,  This  is  a  workable  procedure 
because  if  a  signal  is  present  for  a  relatively  short  time  compared 
to  the  time  taken  for  the  noise  measurement,  then  the  short-time  rise 
in  power  due  to  the  signal  has  a  negligible  effect  on  the  long-time 
a/erage.  If  the  noise  level  is  estimated  in  thiJ  ’•my,  then  by  the 
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central  limit  theorem  the  probability  density  p(N^)  does,  in  fact, 
approach  the  Gaussian,  i.e. 


i  ,  N.  -  N 
1/1  o 


p(Ni) - ± -  e 


N 


/ 2tt 


(70) 


N 


where  N  is  the  mean  value  of  N.  ,  and  oC  hhe  variance  of  the  measure- 

o  l  fi 


ment.  Then 


“  '  2  ’  2  0  (% 


(71) 


where 


Z_  * 
a 


8u  u) 

—a  -a  -  N 

Tu)2  0)^  o 


^2  a 


(72) 


N 


As  before,  if  a  is  specified, 
held  uq  is  given  by 


Tr- 

Ct 


K_, 

a 


a  constant. 


Hence,  the  thres- 


u 

o 


8 


(7  3) 


The  conditional  miss  probability  8(N^)  is  similarly  approximated 
by  a  unit  step  function. 

B(Nl)  -  u(Nl  -  Nu)  (74) 


where  N^  is  the  value  of  noise  level  for  which  -  0.  This  can  be 
obtained  by  replacing  N  by  N. .  in  Eqs.  (49)  and  (50),  and  setting  Z  -  0. 

O  1. 1.  p 

The  value  of  u  to  be  used  in  Eq.  (49)  it  that  of  Eq.  (73)  (with  N  left 
0  o 

unchanged,  since  uq  is  a  fixed  threshold).  The  result  is 


Note  that  rrjj  generally  decreases  with  observation  tin.  ,  the  de¬ 
crease  is,  however,  limited  by  the  degree  to  which  the  sea  noise  is  non- 
stationary . 
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N  +  K-v/2  0„ 
o  &  N 


2 

a“ 

o 


(75) 


Then,  if  p(N^)  as  given  by  Eq.  (70)  is  substituted  in  Eq.  (65),  the 
average  detection  probability  is  found  to  be 


p(D)  =  (1*0) 


(76) 


The  detection  range  h-j  can  be  defined  here  ,»s  the  value  of  h  for 
which  p (D)  *  this  gives 


a)  /  S 

-  k  —  /— L- 

n*  VA.  % 


(77) 


A  curve  of  detection  range  versus  for  the  example  20-knot  target 
used  in  Section  V  is  given  in  Fig.  13;  the  detection  probability  is 
shown  in  Fig.  14.  Note  that  in  the  present  case  both  h-j  and  p(D) 
are  functions  of  K-  .  However,  since  even  for  rather  large  changes 

vX 

in  faise-alarm  rate  K—  only  changes  a  small  amount,  the  effect  of 
change  of  false-alarta  rate,  although  no  longer  completely  negligible 
as  in  Section  V,  still  has  only  a  very  small  effect  on  h^  or  p(D). 
This  is  clearly  shown  in  Figs.  13  and  14. 

In  comparing  the  detection  range  bj  of  this  section  with  h^ 
of  the  last  section,  one  must  note  that  2AN  of  Section  V  is  the 
maximum  uncertainty  range,  while  of  this  section  is  the  standard 
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Fig.  13— Detection  range  versus  aN 


o.oi 
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detection  probability  versus  range 
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deviation  of  the  noise  uncertainty.  For  a  rectangular  distribution, 
the  standard  deviation  would  be 


a 


N 


(78) 


If  this  is  used  in  Eq.  (77),  the  comparison  with  the  detection  range 
of  Section  V  (Eq.  (6G;)  cives 


Thus , 


(79) 


The  detection  probability  drops  off  much  more  rapidly  for  ranges 

in  excess  of  the  detection  range  h—  here  than  in  the  last  section. 

This  is,  of  course,  due  to  the  fact  that  if  the  noise-level  uncertainty 

has  a  Gaussian  distribution,  it  is  much  less  likely  that  the  actual 

noise  level  differs  from  the  mean  value  than  if  the  uncertainty  has 

a  uniform  distribution,  assuming  that  the  variances  are  comparable. 

The  effect  of  transmission- loss  curves  other  than  the  6  db  per 

distance  doubled  curve  considered  here  is  essentially  the  same  as  in 

Section  V.  The  argument  of  ©  in  the  equation  for  P(D)(Eq.  (76)) 

-2 

would  be  proportional  to  g(h/k)  rather  than  h  .  However,  because 
of  the  nonlinear  distortion  of  the  argument  by  the  error  function, 
a  g(h/k)  such  as  that  shown  in  Fig.  12  would  result  in  a  much  more 
violent  fluctuation  of  P(D) ;  in  fact,  it  can  be  concluded  from 
Fig.  14  that  one  would  have  approximately 


Fig.  15 — Average  detection  probability  versus  range 
for  arbitrary  transmission-loss  function  g(h/k) 
and  Gaussian  noise-level  uncertainty 
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VXI .  CONCLUSION 

Several  aspects  of  the  problem  of  detecting  a  submarine  by  a 
single  hydrophone  sensor  have  been  considered.  The  signal  emitted 
by  the  submarine  is  assumed  to  consist  of  broadband  noise  only,  and 
under  this  condition,  it  is  clear  that  the  only  distinguishing  feature 
the  detector  can  use  to  decide  whether  a  target  is  present  or  not  is 
the  increase  in  noise  power  caused  by  the  presence  of  the  targe.. 

A  simple  power  detector,  in  which  the  received  power  is  integrated 
over  the  time  that  the  target  is  within  range  is  therefore  essentially 
equivalent  to  the  optimum  detector. 

In  general,  the  optimum  observation  time  for  a  submarine  moving 
past  the  receiving  sonobuoy  is  inversely  proportional  to  the  velocity 
of  the  submarine.  Hence,  if  the  noise  emitted  by  the  submarine  were 
not  a  function  of  the  velocity,  for  example,  if  it  consisted  exclusively 
of  internally  ge^'rated  machinery  noise,  a  slow  submarine  would  be  some¬ 
what  more  detectable  than  a  fast  one.  However,  even  under  the 
unreasonable  assumption  that  the  background  noise  level  is  known 
precisely,  the  detectability  index  is  proportional  only  to  the 
square  root  of  the  observation  time,  while  for  the  more  realistic 
case  of  unknown  background  noise  level,  the  detectability  index  is 
practically  independent  of  the  observation  time.  ’Thus,  the  increase 
in  detectability  resulting  from  the  increased  time  of  observation  of 
the  slower  submarine  is  extremely  small.  On  the  other  hand,  the 
detectability  index  is  directly  proportional  to  the  power  level  of 
the  signal  whether  the  noise  level  is  known  precisely  or  not. 

Thus,  if  the  noise  emitted  by  the  submarine  increases  at  all 
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with  velocity,  the  increase  in  signal  power  will,  in  practically  all 
cases,  result  in  the  faster  (and  noisier)  submarine  being  more 
detectable.  This  is  especially  true  when  the  background  noise  level 
is  unknown,  in  which  case  this  statement  is  true  without  any 
qualification. 

The  major  factor  limiting  che  maximum  detectable  range  is  the 
uncertainty  about  the  level  of  the  background  noise  level.  If  this 
uncertainty  is  such  that  the  noise  level  is  certain  to  be  somewhere 
between  an  upper  and  a  lower  bound,  then  a  signal  is  in  effect  de- 
tectable  only  if  the  received  power  is  larger  than  the  maximum 
conceivable  background  noise  level.  In  this  case,  one  can  define  a 
minimax  detection  range  such  that  targets  at  a  smaller  range  are 
essentially  perfectly  detectable,  while  the  probability  of  detection 
decreases  sharply  for  ranges  above  this  value.  If  the  transmission 
loss  is  proportional  to  the  square  of  the  distance,  then  the  mlnimax 
detection  range  is  inversely  proportional  to  the  square  root  of  the 
relative  noise-level  uncertainty.  It  is  clear  that  the  performance 
of  the  buoy  depends  critically  on  the  accurate  determination  of  the 
noise  level. 

If  it  can  be  assumed  that  the  noise  level  is  constant,  it  can 
be  estimated  within  the  senobuoy  by  making  a  long-term  power  measure¬ 
ment  of  the  signal  received  by  the  buoy.  If  the  time  of  this 
measurement  is  long  compared  to  the  time  that  target  is  within  range, 
the  overage  power  reading  obtained  is  essentially  that  of  the  noise 
only.  This  method  of  estimating  the  noise  level  results  in  an  un¬ 
certainty  that  can  be  approximately  described  by  a  Gaussian  distri¬ 
bution  whose  mean  is  the  nominal  value  of  the  noise  level.  There 
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are  then  no  upper  or  lower  limits  on  the  possible  noise  level, 
but  detection  still  depends  on  the  received  signal  power  being 
substantially  greater  than  one  would  expect  trom  the  noise  estimate. 
Quantitatively,  if  the  transmission  loss  is  inversely  proportional 
to  the  square  of  the  distance,  the  detection  range  is  proportional 
to  the  square  root  of  the  signal  power  divided  by  the  standard  devi¬ 
ation  of  the  noise  level. 

The  suggested  method  for  internally  measuring  the  noise  level 
results  in  a  two-channel  detector.  One  channel  incorporates  a  long¬ 
term  integrator  having  a  time  constant  Tq  to  estimate  the  noise. 

The  other  one  contains  a  short-term  integrator  with  a  time  constant 
T^  “  3h/v,  the  time  that  the  target  is  within  range.  The  outputs  of 
the  two  integrators  are  compared  in  order  to  arrive  at  a  decision 
as  to  whether  a  target  is  present.  This  system  will  generally  report 
any  rise  in  received  power  lasting  for  a  time  on  the  order  of  as 
a  target;  it  will  tend  to  ignore  changes  in  received  power  that  are 
either  very  much  longer  or  very  much  shorter  than  T^  (unless  the 
increased  signal  powe^  was  quite  large).  The  system  might,  therefore, 
tend  to  reduce  false  alarms  that  are  due  to  anomalous  transmission 
from  very  distant  targets.  On  the  other  hand,  short-term  changes  in 
the  background  noise  level  would  result  in  false  alarms. 

It  is  difficult  to  see  how  any  system  whose  only  basis  for 
detectior  is  the  short-term  change  in  noise  level  caused  by  a  target 
could  be  made  insensitive  to  the  same  sort  of  change  by  the  back¬ 
ground  noise  and  still  be  able  to  detect  targets. 
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The  proposed  two-channel  detector  appears  to  have  certain 
desirable  properties,  and  qualitatively  it  probably  functions  as 
described.  However,  it  is  planned  to  investigate  its  performance 
quantitatively  and  in  more  detail. 
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Appendix  A 

THE  LIKELIHOOD -RAT 10  DETECTOR  FOR  DETECTION  OF  A 
GAUSSIAN  SIGNAL  IN  GAUSSIAN  NOISE  BACKGROUND 

Suppose  that  the  received  signal  has  the  form 

x(t>  -  s(t)  +  n(t)  (A- 1) 

where  s(t)  is  the  signal  that  would  be  observed  if  there  were  no 
noise,  and  n(t)  is  the  noise.  Both  s(t)  and  n(t)  are  assumed  to  be 
Gaus«lan  random  processes  with  zero  mean;  hence  x(t)  is  also  Gaussian 
with  zero  mean.  It  is  always  possible  to  represent  such  a  signal  by 
a  set  of  samples 

X'  -  (A- 2) 

where  the  prime  denotes  the  transpose.  The  kind  of  samples  that  are 
used  depends  on  the  application.  In  some  cases  time  samples  are 
convenient,  i.e.,  x^  ■  x(tp,  ■  x^),  etc.  In  other  cases  it  is 
more  convenient  to  expand  x(t)  in  a  Fourier  series  over  the  observation 
interval  T.  In  that  case  the  elements  of  X  can  be  considered  to  be 
the  Fourier  coefficients.  In  either  case  the  elements  of  X  are  Gaussian 
random  variables  with  zero  mean. 

In  view  of  the  equivalence  between  the  samples  and  the  continuous 
function,  one  can  say  that  the  probability  of  a  particular  realization 
of  x(t)  is  the  joint  probability  that  the  set  of  samples  acquires  the 
particular  value  yielding  this  realization.  Thus,  suppose  that  the 
received  signal  consists  of  noise  only.  The  covariance  matrix  of  the 
elements  of  the  sample  vector  X  is  defined  by 


N-a 


(A-3) 
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where  the  symbol  represents  the  statistical  average  of  the  quantity 

in  the  bracket  conditional  on  the  hypothesis  indicated  by  the  subscript, 
and  where  Q  *  is  assumed  to  exist.  Then  the  probability  density  of  x(t), 
given  that  signal  is  absent,  is  the  probability  density  of  X  given  that 
signal  is  absent,  which  is 


Po<S>  - 


_ 1 _ 

(2it)  n/2  [det  £]1/2 


The  covariance  matrix  of  X  for  x(t)  ■  s(t)  is  defined  by 


(A-4) 


(n')s*i  (4-5) 

Since  signal  and  noise  are  independent,  the  covariance  matrix  for  X  if 
signal  and  noise  are  both  present  is 


a  2'>s+n  -  i  -  s 


Then  the  probability  density  of  x(t)  given  that  signal  and  noise  are 
both  present  is 


- 


(2tt) n/2  [det  (P  +  £)]1/2 


exp|-  “  X^PK^xj-  (A- 6) 


The  likelihood  ratio  is  defined  as  the  ratio  of  p^(X)  to  pQ(X).  It 
is  therefore  given  by 


L(X) 


Pj(2> 

p0<2> 


[ 


det  j) 


det  (P  +  Q) 


1/2 


exp{-  \  X'[(P  +  S)"1-  fi-1]x}  (A-7) 
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In  practice  it  is  generally  reasonable  to  assume  that  the  input 
«ignal-to-noise  ratio  is  very  small,  so  that  the  elements  of  the  P 
matrix  are  all  very  much  smaller  than  the  corresponding  elements  of 
the  £  matrix.  Then  the  ratio  of  the  determinants  in  Eq.  (A-7)  is 
approximately  unity.  Also,  the  exponent  may  be  expanded  as  follows 

-  \  x'[(£  +  a)'1  -  sl1]l  -  - 1  +  mV1  - 1]  x 

-  -  \  x'<f 1  [i  -  P(f 1  +  as'1)2  - 

...  -  i]  X  (A-8) 

"  \  x'SLmlnh.  -  2  x'sfVVV-  •  •  • 


where  £  is  the  unit  matrix. 

If  the  signal  and  noise  are  "white  noise,"  (i.e. ,  with  uniform 
power  spectrum)  then  P  ■  S£  and  Q  **  .  where  S  and  N  are  the  signal 

and  noise  power  respectively.  Then  £  V-Kl)11  s'W1 
-  (  -  ^  I;  etc.  Hence,  the  magnitude  of  the  higher-order  terms  in 
gq.  (A-8)  decreases  with  S/N,  and  for  small  S/N  only  the  first  term 
needs  to  be  considered.  If  X  is  not  "white  noise,"  the  higher- 
order  terms  are  still  negligible,  but  the  demonstration  of  this  fact 
is  somewhat  more  difficult. 

Generally,  it  is  more  convenient  to  deal  with  the  logarithm  of 


the  likelihood  ratio  than  with  the  ratio  itself.  This  causes  no 
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difficuHy  ",.lnce  the  logarithm  is  a  monotonic  function  so  that  if  the 
likelihood  ratio  exceeds  a  given  threshold,  its  logarithm  exceeds  a 
different  threshold.  For  small  signal-to-noise  ratio,  the  above  dis¬ 
cussion  indicate  s  that 

log  L(X)  *  U  -  j  x'if  *lSf  **  (A-9) 

The  quantity  U  is  the  test  statistic  of  the  problem,  and  it  is 
compared  to  the  threshold  U  to  decide  whether  or  not  a  signal  is 
present.  It  is  essentially  the  quantity  that  must  be  computed  by 
the  optimum  receiver.  Although  the  likelihood  ratio  is  a  ratio 
of  probability  densities,  U  is  obtained  from  X  by  a  deterministic 
transformation  of  the  vector  X  representing  the  received  signal.  U  is  a 
random  variable  _ecause  the  vector  X  is  a  random  variable. 

If  U  >  Uq  the  decision  is  that  a  signal  is  present.  The  con¬ 
ditional  false-alarm  probability  a  is  therefore  the  probability  that 
the  randem  variable  U  exceeds  the  threshold  UQ,  given  that  the  signal 
is  actually  absent,  i.e. 

OD 

a  -  f  p  (U)  dU  (A- 10) 

JU  0 
o 

where  PQ(U)  is  the  conditional  probability  density  of  U  when 
there  is  no  signal.  Similarly,  the  conditional  probability  of  a 
correct  detection,  (1-3),  is  given  by 

f® 

(1-3)  -  P,(U)  dU 

u  1 
o 


(A-ll) 
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where  p^(’J)  is  the  probability  density  of  U  when  there  is  a  signal  present. 

In  order  to  evaluate  a  and  (1-g) ,  Pp(U)  -nd  p^(U)  must  be  computed. 

This  is  in  general  rather  difficult  because  the  transformation  given  in 

Eq.  (A-?}  is  nonlinear.  If  the  matrix  Q  ^PQ  ^  equalled  DI,  where  D  is 

a  constant  and  1^  is  the  unit  matrix,  U  would  reduce  to  the  sum  of  squares 

of  the  elements  of  X.  Then,  since  X  is  a  Gaussian  vector,  the  probability 

•k 

density  of  U  would  be  a  chi-square  density  in  n  «  2TW  degrees  of  freedom. 

-1-1  ** 

For  general  £  PQ  ,  the  probability  density  of  U  is  what  Middleton 

refers  to  as  a  "generalized  chi-square"  density.  It  can,  however,  be 

shown  (Ref.  1,  Section  17.2-1,  and  Ref.  8,  Section  20.2)  that  if  n  is 

large,  and  if  the  eigenvalues  of  the  ma;rix  ^  ^PQ  *  are  small  (which 

means,  in  effect,  that  the  signal- to-noise  ratio  is  small),  both  PQ(U) 

and  p^(U)  approach  the  Gaussian  form.  Both  of  these  conditions  would 

normally  be  expected  to  hold  very  well  for  problems  of  interest  here. 

It  therefore  will  be  assumed  that  PQ(U)  and  p^U)  are,  in  fact,  Gaussian. 

A  Gaussian  distribution  is  completely  determined  by  its  mean  and 

variance.  The  mean  of  PQ(U)  is  given  by 


n  n 


2  \  5  X  "  2  I  I  aiJ  (  Vj  \ 

1-1  J-l 


(A-12) 


where  a^  is  the  general  element  of  (}  *~PQ  However, 


(vj)n 


(9) 


ij  * 


the  general  element  of  (J>  by  definition.  Also,  it  can  be  shown 

that  the  double  sum  £  Zb,,  c.,  of  the  product  of  all  the  elements 

i  j  J  J 

of  the  two  matrices  B  and  C  is  equal  to  the  sum  of  the  diagonal  terms, 
or  ti.ee,  of  the  product  matrix  BC,  written  tr (BC) .  Hence 


See  Ref.  8,  Section  18.1. 

Ir* 

See  Ref.  1,  Section  17.2-1, 
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u  u 


n  i 

i-i  j-i 


L  aij  qij 


~  trCs"1^"1^) 


(A- 13) 


-  2  trC <f  l£)  -  J  tr(2£l) 


Similarly,  one  finds  that  the  variance  of  U,  given  that  noise  only 


is  present 


I (ci'a'W  )n  -  ;( xW1* )k 


n  n  n  n 


-  4  l  I  I  l  \i\l  {  VjVj  )  -  *  ]  <A-“> 

i-1  j-1  k-1  J&-1 

where  the  a^  and  a^  are  again  general  elements  of  £  Since 

the  x,  are  Gaussian  variables,  one  can  expand  the  quadruple  average 

X 

into  products  of  covariances. Therefore 


n  n  n  n 


l  1  I  l\s\t(  Vj Vi  \ 

i-1  J-i  k-i  *■! 


N 


n  n  n  n 


III  l  ‘i)\t  [  (  Vj  }N  (  Vi  >N 

t _ 1  l _ 1  a _ 1  •’  Al 


i-1  J-1  k-1  f-1 
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Thus 


n  n  n  u 


l  1  I  l  “u  \t  ( VjYi  > 

-1  4-1  t-1  #«1  N 


i-1  j-1  k-1  je-1 

n  n 


[  x  x 


i*l  j-1 


+  2 


n  n  n  n 

III  ZatJaw 

i-1  j-1  k-1  4-1 


Xi*k 


[ 


tr 


+ 


L. 


Finally 


(A- 15) 


(A- 16) 


Similarly,  it  can  be  shown  that  the  mean  ^  of  p^U)  is  given  by 
“  2  (  LP2  LX  ^  J  tr[(P2-1)  ]  +  \  tr(P(f  *)  (A-17) 


and 


2  1 
al  '  2  tr 


-1  2  -1  2 

m  >  a  +  £2  > 


\  tr[(SS"1)  J  - 


2.. 


(A- 18) 


where  the  approximation  implies,  as  before,  that  the  signal-to-noise 
ratio  is  small  so  that  the  elements  of  P£  1  are  all  very  much  less 
than  unity. 
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From  Eqs.  (A-10),  (A- 13),  and  (A-16) ,  the  false-alarm  probability 


is  found  to  be 


1  1  « 
a  «  2  "  2  ® 


Uo  -  \  trCP^'1) 


(A- 19) 


where 


o  z  _t2 

0  (z)  =  J  e  dt 


Similarly,  from  Eqs.  (A-ll),  (A-17),  and  (A-18),  the  probability  of 


a  correct  detection  is 


1  !  .  ®.  ■  1  “OS'1) 

d-B)  -  2  -  2  ®  ^ . — 


r[£2"1)2 


(A- 20) 


Comparison  of  the  arguments  of  the  ®  functions  in  the  two  cases 
indicates  that  the  difference  is  the  quantity  d//2,  where 


/ l  [  -i  21 

2  trt^  }  - 


(A-21) 


d  is  referred  to  as  the  "detection  index"  and  may  be  considered  as  the 
figure  of  merit  of  the  detection  system.  For  a  given  false-alarm 


probability  which  fixes 


[Uo  ’  2  tr<£9  1)]y^trC(P2  l)2]>  the  largei 


d  is,  the  greater  is  the  probability  of  true  detection.  Inspection  of 
Eqs.  (A- 13),  (A- 16),  and  (A-17)  indicates  that 


»*i  " 


(A-22) 


The  optimum  receiver  forms  U  from  x(t)  according  to  Eq.  (A-9). 
One  can  think  of  pj  and  as  the  "dc"  component  of  the  output  under 
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the  hypothesis  that  the  signal  is  present  or  absent  respectively,  and 

a  rv  c,  is  the  ms  fluctuation  of  U.  Then  d  is  the  ratio  of  the 
o  1 

difference  of  useful  output  from  the  detector  to  the  output  fluctu¬ 
ation,  i.e.,  it  is  like  an  output  signal-to-noise  ratio.  This 
suggests  that  this  same  figure  of  merit  be  used  in  evaluating  sub¬ 
optimum  systems. 
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Appendix  B 

DETECTION  WHEN  SIGNAL  OR  NOISE  POWER  IS  UNKNOWN 


In  the  previous  discussion  it  has  been  tacitly  assumed  thac  the 
signal  power  and  noise  power  are  known;  otherwise  the  matrices  P  and 
(J,  which  are  proportional  to  signal  and  noise  power  respectively, 
could  not  be  completely  determined.  In  practice  neither  of  these 
power  levels  would  be  known  exactly. 

Consider  first  the  problem  of  unknown  signal  power.  It  seems 

clear  on  heuristic  grounds  that  the  false-alarm  probability  a  should 

not  be  a  function  of  the  signal  power,  and  this  is  corroborated  by 

■* 

statements  in  the  literature.  However,  it  appears  from  Eq.  (A-19) 
that  a  does  depend  on  both  P  and  £  and,  therefore,  on  both  signal 
and  noise  power.  In  order  to  resolve  this  apparent  discrepancy,  we 
note  that  the  test  statistic  U  defined  in  Eq.  (A-9)  is  generated  from 
the  received  signal  by  means  of  a  processor  which  in  some  way  realizes 
the  matrix  operation  ^PQ  ^ .  Since  the  elements  of  P  are  proportional 
to  the  signal  power  S,  the  "gain"  of  the  processor  is  proportional  to 
S.  This  is  the  only  reason  why  U  is  proportional  to  S;  the  received 
signal  x(t)  (consisting  of  noise  only)  has  nothing  to  do  with  it.  The 
threshold  Uq  for  a  desired  false-alarm  probability  therefore  aJso 
depends  on  S  only  because  the  gain  of  the  processor  does.  If  the 
processor  gain  is  arbitrarily  changed,  the  false-alarm  probability  q 
will  remain  unchanged  if  the  threshold  Uq  is  changed  in  the  same 
proportion  as  the  processor  gain.  Thus,  it  is  possible  to  base  both 
the  processor  gain  and  the  threshold  on  a  normalized  £  matrix,  i.e., 

*See  Ref.  2,  pp.  133-134. 
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to  set 


l  ‘  s  £ 

ad  (B-l) 

U  «  S'J' 

o  o 

Substitution  of  Eq.  (B-l)  into  Eq.  (A- 19)  then  indicates  that  a  is 
independent  of  S. 

This  same  argument  does  not  hold  for  unknown  noise  level.  In 
order  to  show  this,  suppose  that  the  matrix  £  is  also  normalized,  i.e. 

£  -  N  3  (B-2) 


and  that  a  test  statistic  u  is  formed  by  a  normalized  processor 
5  1£3  *  which  is  now  independent  of  both  signal  and  noise  power,  i.e. 

u  -  Ix's'W1*  (B-3) 


It  is  easily  seen  that  for  x(t)  consisting  of  noise  only,  both  the 
mean  and  standard  deviation  of  u  are  proportional  to  N  through  the 
received  signal  vector  X.  If  u  is  compared  to  the  threshold  uq, 
the  false-alarm  probability  is  given  by 


a 


1 

2 


uc  "  2  N  tr(pqX) 
N  N/tr[ (pq  V] 


(B-4) 


This  is  clearly  dependent  on  N.  The  noise  level  canoot  be  normalized 
out  because  it  depends  on  the  received  signal  rather  than  being  simply 
a  gain  adjustment  in  che  processor  as  is  the  signal  level. 
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Suppose  now  that  the  received  signal  x(t)  consists  of  noise 
with  a  known  power  level  N,  and  signal  with  an  unknown  power  level 
S,  and  that  this  signal  is  processed  through  a  normalized  processor 
as  in  Eq.  (B-3).  For  small  signal- to -noise  ratio  (which  can  be 
assumed  to  exist  even  though  the  exact  signal  power  is  not  known), 
the  variance  of  u  is  the  same  as  with  noise  only  present.  However, 
the  change  in  the  mean  due  to  the  presence  of  the  signal  is  now 


(  “L  '  (  U  )„*  2  s  tr(Pq'1) 


(B-5) 


S+N  '  'h 

Therefore,  the  detection  index  d  now  takes  the  form 


N  V  2 


tr(Ea”1) 


(B-6) 


and  Eq.  (A-20)  for  the  probability  of  correct  detection  is  changed 
to 

(1-8*  '  ^  “  8 

where 

u 


is  the  argusent  the  error  function  defining  a. 

It  should  be  noted  that  although  in  this  discussion  the  matrices 
£  and  £  were  normalized  with  respect  to  total  signal  and  noise  powers, 
any  other  normalization  proportional  to  signal  and  noise  power  would 
yield  identical  results. 


k(or) 


-i5  / 
2  y  v 


tr(pq  ) 


(B**7) 


-  \  N  tr(23  l) 


7 


(B-8) 


i  ytr(£a  l) 
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